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Abstract 

We present an action principle formulation for the study of motion of an 
extended body in General Relativity in the limit of weak gravitational field. 
This gives the classical equations of motion for multipole moments of arbitrary 
order coupling to the gravitational field. In particular, a new force due to 
the octupole moment is obtained. The action also yields the gravitationally 
induced phase shifts in quantum interference experiments due to the coupling 
of all multipole moments. 
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The study of motion of extended bodies possessing multipole moments in the gravita- 
tional field has a long history [1]. The starting point was the Einstein - Infeld - Hoffman 
[2] derivation of the geodesic equation for a point test particle from the gravitational field 
equation and the conservation law for stress energy tensor. The test particle approxima- 
tion breaks down if body's extension in space is non negligible compared to the radius of 
curvature of the background field and secondly when the back reaction due to the body on 
the background field is non ignorable. In this letter, we shall be concerned with the former 
aspect. This is particularly motivated by the fact that astrophysical bodies like planets and 
stars are extended and should in a realistic analysis be treated as such. The interaction of 
covariant generalization of Newtonian multipole moments with the gravitational field will 
be given by their coupling to Riemann curvature and its derivatives. This would appear as 
modification to the geodesic equation. 

The modification to the geodesic for a spinning body is given by the Mathisson- 
Papapetrou equation [3,4], which may be extended to a particle with intrinsic spin [5]. 
Subsequent to the treatment of spinning bodies, various authors have obtained the correc- 
tions up to the covariant generalization of Newtonian quadrupole moment [6-8]. A com- 
prehensive study of the problem including comparison of various approaches and results is 
carried out by Dixon [8], but, to our knowledge, no one has obtained corrections to geodesic 
equation arising due to coupling of covariantized higher order Newtonian multipole moments 
with gravitational field. More importantly a procedure to derive equations of motion of ex- 
tended bodies, with arbitrary multipole moments, through an action principle has not been 
obtained during the past 65 years in which the equations of motion in a gravitational field 
have been studied [9]. In the absence of a general principle to obtain the action made up of 
terms uniquely attributed to couplings of all multipole moments, such a task is very difficult. 
This is precisely what we wish to do in the following. 

While the action is used in classical physics only as a tool to obtain the equations of 
motion, the action is directly observable in quantum physics as the phase of the wavefunc- 
tion. Therefore, the phase shift produced by the coupling of multipole moments with the 
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gravitational field can in principle be measured (action giving an algorithm to calculate 
various multipole phase shifts). Neutron interferometry provided the first instance where 
effect of the earth's gravitational field on the phase of the neutron wavefunction was ob- 
served [10]. Interesting gravitational analogs [11] of the topological Aharanov-Bohm [12] 
and Aharanov-Casher [13,14] phase shifts have also been proposed. 

In this paper we present a formalism which yields in a simple and elegant way the 
corrections to the geodesic equation up to all multipoles of the extended body. The equations 
of motion for multipoles simply follow from variation of our action. As a demonstration of our 
formalism we obtain for the first time corrections to the geodesic equation till the covariant 
generalization of Newtonian octupole moment. Moreover, our action gives the quantum 
phase shifts in interferometry due to the coupling of all multiple moments to the gravitational 
field. Our formulation of action principle for extended body may facilitate and prompt 
further investigations in the important and emerging area of interface between quantum and 
gravitational realms [15]. Particularly, the ongoing experiments in atomic [16], molecular [17] 
and Bose-Einstein condensate interferometry [18] hold promise for experimentally testing the 
new gravitational phase shifts that will be obtained in the present letter. 

We envision an extended rigid body as a thin world tube in spacetime and its thickness 
is small compared to the scale over which curvature varies. We would further assume that 
there are no external or internal forces acting on this body apart from the gravitational field 
in which it is propagating. In the thin world tube we choose a reference world-line (z^) 
having 4- velocity = dz^/ds = (1,0,0,0) and define multipole moments with respect to it 
on a spacelike hypersurface. The multipole moments of order 2n are defined as [8] 



where 5%^ = — z^, T^ u is the energy- momentum tensor, and the integration is over the 
spacelike hypersurface identified by the unit normal vector field w a . The above multipole mo- 
ments are defined in a class of coordinate systems that are related by linear transformations 
in order that the expression (1) is covariant. But once they are defined this way, t K1 - KntMU can 




(1) 
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be transformed to any arbitrary coordinate system as a tensor. All the relativistic equations 
in the present letter are covariant with respect to the above linear transformations, if not 
with respect to general coordinate transformations. 

Since what we do must be consistent with the Newtonian theory in the appropriate 
limit, we shall now establish the relation between the covariant multipole moments of order 
2n and the covariant generalizations of the anti-symmetric spin tensor , the symmetric 
quadrupole moment tensor J* 3 ' and the symmetric octupole moment tensor 0* jfc in Newtonian 
gravity. We Write the Newtonian potential energy U of the body with mass density p(x) in 
terms of the potential (f>(x) expanded in a Taylor series around the central world-line, 

U = J p(x)<f ) (x)d 3 x = m<j ) (z) + d l d l( f ) \ z + ^P'd l d j( j ) \ z + ^O^ k d z d J d k( j ) \ z + ... (2) 

where the mass m — J p{x) d 3 x, the dipole moment d l — J p{x) 5x l d 3 x, the quadrupole 
moment P^ = J p(x) 5x l Sx^d 3 x and the octupole moment 0^ k = J p(x) 5x l 5x^ 5x k d 3 x, with 
Sx i = x i - z\ In view of eq.(l), we identify f 00 = d\ = , t ijk0 ° = O ijk . The Spin 

tensor(orbital angular momentum) in the Newtonian limit is defined as S lj = 2 / pSx^v^cPx 
where v l = dx % jdt. The Spin tensor then satisfies dS^/dt = 2p^u^ — 2/ p5x^d^(pd 3 x, where 
u % = dz % jdt and the momentum p l — J p(x) v l d 3 x. Using Taylor expansion of the potential, 
and choosing z l to be the center of mass so that d % = 0, the spin propagation equation up 
to octupole term becomes 

^-S ij = 2pM - 2 I k ^d k <f)\ z - O kr[l d l] d k d r ct)\ z (3) 
dt 

The covariantization of this spin tensor leads to 

In the weak field limit, the metric is g^ v = r]^ v + h^ v with h^ u <C 1. In the Newtonian limit 
g 00 = 1 + 2(f) which implies = h 00 /2. Thus the dipole term in eq.(2) leads to the covariant 
form t MQ/3 h af 3^\ z , where t^ a/3 includes the spin tensor. We shall now choose the reference 
world line to be the center of mass so that d l = 0. We consider only matter distributions 
for which [19] 



ir a = s^ u u a \ (5) 

whereS^ satisfies S^u v = 0. 

The quadrupole term leads to didj(f)\ z I lj = ^didjh 00 \ z P j = — Roioj\ z We thus have 
di djhoo\ z t 1 - 700 = — 2 Roioj\ z I 1 ^ which covariantizes to 

Kp^\ z t^ = -2R a ^ v \ z I^u a u p (6) 

where the covariant quadrupole tensor, = I up . Similarly, the octupole term leads to 

K^ v Xt^ = -2 R a ^ a \ z u a v? (7) 

where 0^ ua is the fully symmetric covariant octupole tensor. Eqs.(5), (6) and (7) are the key 
relations which would unambiguously provide the connection between covariant multipole 
moments with their Newtonian analogs. 

In the Newtonian approximation the phase shift in quantum mechanical interference due 
to the gravitational field is 

l r „ . l , r , . /■„■„.. r l 



*= s 



Judt = --(Jm</>dt + J cf d^dt + ^ - / ij <9; 9,- 0dt + J- O ijk di dj d k <pdt + ...) 

(8) 



where U is given by the expansion in using eq.(2). The first term of eq.(8) corresponds 
to the phase shift observed in the COW experiment [10], and the subsequent dipole and 
quadrupole terms are corrections to it. The higher order multipole contributions to the 
phase shift may also be obtained from this expansion. In General Relativity this phase shift 
is obtained by letting the path ordered operator resulting from the covariant generalization 
of this Newtonian phase shift act on the initial wavefunction. Using eqs.(2, 5, 6, 7) and 
noting that in the linear field limit uj a b pS b a = h a/3 ^S ap , where uj a b p are Ricci rotation 
coefficients,this path ordered operator is given by 

g = Vexp [ - % - J (- m + l - uj\pS a b v?- X - R a ^J a/3 u»u v - l - J R a ^ v , p O a ^ u" u v + ...)ds 

(9) 



In the special case of I af3 and tjk being zero, this result is in agreement with the gravita- 
tional phase shift for intrinsic spin [5]. 

We shall now obtain this expression in the weak field limit of General Relativity starting 
from the action principle. Choose a coordinate system such that along the reference world- 
line g^ v = rj^y. We would here like to recall that the extended body under consideration 
is rigid and subject to only external gravitational force and no other forces (external or 
internal). In that case the monopole moment corresponds to the mass (m) of the body with 
whole matter concentrated on the reference world-line and ignoring the back reaction on the 
background gravitational field we write m = p a u a . The higher order multipole moments 
are defined by considering the metric perturbations as one moves away from the reference 
world-line z^. For simplicity we restrict perturbations of the metric to the first order [20], 
and later covariantize the equation of motion by changing ordinary derivative to covariant 
derivative. The action 

S = j V^9 Cd 4 x = j C\ g ^=^ v d 4 x + ^ v/=# Cd 4 x- j £\ g ^= v ^ d 4 x^j . (10) 
The first term is the kinetic energy term, / y/^^g C\ 9lll/=Vflv d 4 x = — J p a u a ds. So, 
S = - J Pa u a ds + J 5g, v ^-(y/=jj£) d 4 x + ... = - j ' Pa u a ds+ l - J 5g, u ^ T» v d 4 x + ... 

(11) 

Note bg^y = g^ u — rj^ = h^ u and hence it can be written as 5g^ u = h^ Uta \ z 5x a + 
\ hn U;a p\ z 5x°~ 5x p + ...where we have h pv (z) = 0. In the present weak field limit, we ne- 
glect all terms that are quadratic or higher order in metric perturbations. Substituting for 
5g^ v in eq.(ll) and using eqs.(l,5,6,7), we finally obtain 

S = - J p a u a ds+ l - J h a ^ S» a u?ds- l -J R a ^ v I al3 u »u»ds- 1 - J R a ^ p O al3p u" u u ds 

(12) 

up to octupole (all derivatives are evaluated on z^). In the linear field limit, since h a/ 3 :fl S m = 
uj a b(3S a b , hence the accumulation of infinitesimal phases arising from eq.(12) is the same as 
what is obtained from the path ordered operator eq.(9) in this limit. 
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We now obtain the equation of motion by extremizing the action (5S = 0) and requiring 
that coordinate variations vanish at end points of the path. This leads to the equation of 
motion 

dp 



ds 



-R^,„ «" r {p « A) + j R^„„ W 0<"><" u» - i ^ (O"' v.") (13) 



In a coordinate system in which h a p^ = 0, (13) gives 



+^ WpRatf* u a 0**> u - I Vp (O^ « a ) . (14) 

Eq.(14) is generally covariant and thus is valid in every coordinate system. Since the dipole 
moment couples to Riemann curvature, it is expected that quadrupole should couple to 
its first derivative and octupole to its second derivative. However there is a coupling of 
quadrupole with Riemann curvature in the second term and a coupling of octupole with 
first derivative of Riemann curvature in the last term which suggests that p a should be 
suitably redefined [21] as, p* = p a — R a \ pu J Am u v — |V ' p R a ^uaO llup u a : , which would then 
yield the expected form 

^ = \ RaX^u U X + l ~ VaRa^u ^/ q][/ V ] + ^ V p V CT R a ^ v U^O^U® . (15) 

This is the equation of motion for a body possessing dipole, quadrupole and octupole mo- 
ments in a gravitational field. The spin propagation equation, eq.(3), can also be covariantly 
generalized to 

E-S ap = 2p^ a - 2 R^^ u^I^ v u a] - R [a ^ a - P u [fl PMv w" 1 (16) 

which on the redefinition of momentum vector form p a to p*modifies as follows: 

_£_ S af3 = 2 p*[* u P] + 2 (R a ^ a u^I^ u u a] - R^ ua u [a I^ [v u a] ) 

-2R^^ a u^I^u^ -' \R [a ^ p u^O^ [u u° ] . (17) 



To simplify the notation we would further define J^ a/3v := -Su^I^u^ and G^ av : = 
-u^O a](3[a u^ and then finally obtain 

Dr>* 1 1 1 

= g Rax,u u x S>" + - J^ u V a R» a p u + - V P V CT R a ^ u G a ^ . (18) 

£_ S af3 = 2p *[« u f3] _ 4 R[a ^ _ 5 ^ [Q ^ G/3]Mp , CT (lg) 

Thus we have obtained for the first time the correction to propagation equations till the 
coupling of covariant generalization of Newtonian octupole moment with the background 
gravitational field. These equations agree with the earlier results obtained till quadrupole 
[8] but the force due to the octupole moment is a new result. Our procedure can be easily 
generalized easily extended to obtain further corrections from due to all higher multipole 
moments. 

This is a very simple and elegant method of deriving the equation of motion for an 
extended body incorporating coupling of multipole moments of arbitrary order with the 
gravitational field. And above all we have found the action for such a body which we have 
also used to compute the gravitational phase shifts in its quantum wavefunction due to the 
multipoles. We hope that the new gravitationally induced phase shifts would be measured 
in future interferometry experiments based on atomic, molecular and Bose-Einstein conden- 
sates [16-18]. Our novel algorithm also yields modifications to the geodesic equation which, 
with present day high precision astronomical observations yielding multipole moments of 
planetary or stellar bodies, can be applied to obtain aberrations in their orbits. It would be 
interesting to relax the other aspect of test body character; i.e. non namely ignorability of 
back reaction, and then study the motion in full generality. 
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